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Abstract — 

We study the phase diagram of the sine circle map lattice with random 
initial conditions and identify the various types of dynamical behaviour 
which appear here. We focus on the regions which show spatio-temporal 
intermittency and characterise the accompanying scaling behaviour. Di- 
rected percolation exponents are seen at some points in the parameter space 
in the neighbourhood of bifurcation boundaries. We discuss this behaviour 
as well as other types of behaviour seen in the parameter space in the con- 
text of the phase diagram. 
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I. Introduction 

THE existence of spatio-temporal intermittency wherein 
laminar regions, which exhibit regular dynamics in space 
and time, co-exist and propagate together with regions which 
show irregular or chaotic bursts, which can be termed turbu- 
lent behaviour, can be seen in a wide variety of spatially ex- 
tended systems. Examples of such systems range from oscil- 
lators, systems which show pattern formation 1 1 1, chemical re- 
actions G), to turbulence @ and fluid flows 0, 0, 0. The 
phenonomenon of temporal intermittency has been studied ex- 
tensively and is relatively well understood |7|. The nature of 
intermittency in spatially extended dynamical systems, however 
has not been understood very well. The presence of spatial as 
well as temporal intermittency has implications for understand- 
ing the physics of pattern formation and for understanding the 
ubiquitous presence of structures in chaotic systems. 

Coupled map lattices, i.e. systems with continuous vari- 
ables which evolve on discrete space-time, are particularly sim- 
ple paradigms for systems with extended spatial dimension 
and show a wide range of interesting dynamical behaviour [ 1 1 . 
Spatio-temporal intermittency has been observed in a variety of 
coupled map lattices |8|, [9|, 111 01 . Ill II . Phase transitions to 
spatio-temporal intermittency are a topic of current interest, and 
the identification of the universality class of this transition has 
led to much discussion in the literature. Specifically, the type 
of spatio-temporal intermittency where there is no spontaneous 
creation of bursts, and a given laminar site can only become tur- 
bulent if infected by turbulent neighbours, has been conjectured 
to lie in the same universality class as directed percolation, with 
the laminar sites being identified as the 'dry' or absorbing states, 
and the turbulent sites being the 'wet' or percolating states. In 
another type of spatio-temporal intermittency, laminar sites can 
spontaneously become turbulent. These two types of intermit- 
tency have been called Type I and Type II spatio-temporal inter- 
mittency in the literature, respectively. 

Recently, the existence of spatio-temporal intermittency has 
been observed at certain parameter values for the sine circle map 
lattice, a popular model for the behaviour of mode-locked oscil- 



lators i fTTl . However, the detailed phase diagram of this model 
has not yet been obtained. In this paper we study the phase 
diagram of this model and identify the various types of dynam- 
ical behaviour which appear here. We next concentrate on the 
regions which show spatio-temporal intermittency and charac- 
terise the accompanying scaling behaviour. Directed percola- 
tion exponents are seen at some points in the parameter space 
in the neighbourhood of bifurcation boundaries. We discuss this 
behaviour as well as other classes of dynamical behaviour seen 
in the context of the bifurcation behaviour as well as the phase 
diagram of the system. 

II. The Model 

In this paper, we consider the coupled circle map lattice 1141 
which has been used to model mode-locking behaviour. The 
coupled sine circle map lattice is defined by the evolution equa- 
tions I1 T2I : 

x n+1 (i) = (l-e)/(a; n (i))+|[/(a; n (i+l))+/(x n (i-l))](niod 1) 

(1) 

where t is the discrete time index, and i is the site index: i = 
1, . . . L, with L being the system size. The parameter e gives 
the strength of the diffusive coupling between site i and its two 
neighbours. The local on-site map is given by 
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where the parameter K is the strength of the nonlinearity, and 
f2 is the frequency of the single sine circle map in the absence 
of nonlinearity. This CML has been studied extensively and 
has a rich phase diagram with many types of attractors and 
strong sensitivity to initial conditions 1 121 . 11 31 . In particular, 
this system also has regimes of spatio-temporal intermittency 
(STI) when evolved in parallel with random initial conditions 
111 II . For weak values of nonlinearity, i.e. for K = 1.0 the 
system shows spatio-temporal intermittency with a unique ab- 
sorbing state at the parameter values SI = 0.064, e = 0.63775, 
and 57 = 0.068, e = 0.73277. These mark the transition from 
a laminar phase to STI. The laminar phase here corresponds 
to the synchronised fixed point x* = ^- sin _1 (^^), and the 
turbulent sites are those which are different from x* . Spatio- 
temporal intermittency has also been observed for high val- 
ues of nonlinearity K — 3.1, this time with infinitely many 
absorbing states. An earlier study of the inhomogeneous lo- 
gistic map lattice had shown spatio-temporal intermittency in 
the neighbour-hood of bifurcations from the synchronised fixed 
point of the system. It is therefore worthwhile to investigate 
the detailed phase-diagram of the present system, identify var- 
ious types of dynamical behaviour, and correlate the observed 
behaviour, especially spatio-temporally intermittent behaviour, 
with the known bifurcations that occur in the system. 
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Fig. 1 

The phase diagram in (e, Q) space for the coupled circle map lattice after 
15000 iterations for a lattice of 1000 sites. Random initial conditions were 
used. The regimes of spatially synchronised and temporally fixed solutions (+) 
and the cluster solutions, where the number of clusters is less than 10 (□) are 
shown . The points at which DP-like exponents are obtained are also marked 

(O). 



III. The phase diagram 

We study the system with random initial conditions at the pa- 
rameter value K = 1.0 as detailed phase diagrams with other 
classes of initial conditions are available at this parameter value 
1121 . 1131 . 1141 . The single circle map shows temporal period 1 
solutions in the region < £1 < and the coupled sine circle 
map shows synchronised, period 1 solutions over the same range 
when evolved with synchronised initial conditions 1 14 1. There 
is a similar region between 1 — < 1. We investigate this 
region for < e < 1 with random initial conditions. The phase 
diagram for this region can be seen in Fig. 1. The following 
types of dynamical behaviour can be found here. 

1. Synchronised solutions: These constitute the dominant be- 
haviour over the bulk of the parameter space and are very stable 
to perturbation. These are indicated by plus signs in Fig. 1 . 

2. Cluster solutions: There are regions which correspond to 
clusters where the variables on the lattice take up a finite num- 
ber of values, but at random sites in the lattice, e.g. a 2— cluster 
is a solution where the lattice variables take up only two values, 
and sites which have the same value for the variable, belong to 
the same cluster. Regions which have less than 10 clusters are 
indicated by boxes in the phase diagram. 

3. Spatio-temporally intermittent solutions can be seen at the 
parameter values marked with diamonds in this region. We dis- 
cuss these solutions in detail in the next section. 

4. Spatio-temporally disordered solutions can be seen in the 
white regions in the parameter space. 

It is clear that the synchronised solutions change to clus- 
ter solutions at some of the phase boundaries, and to spatio- 
temporally disordered regions at others. There are some regions 
where a partially ordered phase is seen in the phase diagram 
where there are many clusters. The boundary between this par- 
tially ordered phase and the spatio-temporally disordered phase 
is not clear. While there are various interesting aspects to the 
behaviour of the cluster solutions, we will concentrate on the 
behaviour of the spatio-temporally intermittent solutions in this 




Fig. 2 

A space-time plot for the spatio - temporal Intermittency seen in the coupled 
circle map lattice at f! = 0.049, e = 0.8495. Horizontal axis is space and time 
is along the vertical axis. The system size chosen is L = 1000 and 15000 
iterations have been discarded. The black region corresponds to the laminar 
state and the white region corresponds to the turbulent state. Every 10 th iterate 
has been plotted. 



paper. 

IV. SPATIO-TEMPORAL INTERMITTENCY AND DIRECTED 
PERCOLATION EXPONENTS 

Spatio-temporally intermittent regions can be seen at several 
points in the parameter space. The points at which spatio- 
temporal intermittency of the first type are seen are marked on 
the phase diagram. As mentioned earlier, £1 = 0.064, e = 
0.63775, and Q = 0.068, e = 0.73277 are two of these 
points with a unique absorbing state with the value x* = 
t^t sin -1 {^jp-). The same kind of spatio-temporally intermit- 
tent behaviour is found at two more points in the parameter 
space viz. il = 0.049, e = 0.8494 and fl = 0.073, e = 0.4664. 
To verify whether the spatio-temporal intermittency here be- 
longs to the universality class of directed percolation, set of 
critical exponents which describe the scaling behaviour of the 
quantities of physical interest have been calculated. The physi- 
cal quantities of interest for such systems are (a) the escape time 
r, which is the number of time steps elapsed before the system 
reaches its laminar state and (b) the order parameter, m(e, L, t), 
which is the fraction of turbulent sites in the lattice at time t. 
From finite-size scaling arguments, it is expected that r depends 
on L such that 

{log L laminar phase 
L z critical phase 
exp L c turbulent phase 

Here, c is a constant of order unity, and the critical point is iden- 
tified as the set of parameter values at which r shows power-law 
behaviour with z being the critical exponent. At e c the critical 
value of the parameter e (other parameters being held fixed), the 
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Fig. 3 

The escape time, r vs the system size, L is plotted at e values below, at and 
above criticality at Q = 0.049 on a log — log scale. The values of e plotted are 
0.845, 0.849, 0.8494, 0.85, 0.852. Power law behaviour is seen at 
e c = 0.8494(+). The data at e c is fitted to an exponent 1.578. 



order parameter m(e,L,t) scales as 

m~{e-e c f, e - e+ (3) 

when the critical line is approached from above . 

The order parameter is expected to obey the scaling relation 

m( ec ,L,t)«i-^ z (4) 

for t « t. The exponent v can thus be extracted once (3 and 
z are obtained from the relations above. 

To extract further critical exponents, we obtain the correla- 
tions from the pair correlation function given by: 

1 L 

Cj(t) = j;22 < ^K+iW > - < x i(t) > 2 (5) 

i=l 

where the brackets denote the averaging over different initial 
conditions. At criticality one expects an algebraic decay of cor- 
relation 1 13): 

where rf is the associated critical exponent. 

The above set of exponents constitute the static exponents of 
the model. The dynamic exponents are given by the spreading 
exponents defined as follows. We consider temporal evolution 
from initial conditions which correspond to an absorbing back- 
ground with a localised disturbance, i.e. a few contiguous sites 
which are different from an absorbing background. The quanti- 
ties of interest are, the time dependence of N(t), the number of 
active sites at time t averaged over all runs, P(t), the survival 
probability, or the fraction of initial conditions which show a 
non-zero number of active sites (or a propagating disturbance) 
at time t and the radius of gyration R 2 (t), which is defined as 
the mean squared deviation of position of the active sites from 
the original sites of the turbulent activity, averaged over the sur- 
viving runs alone. At criticality, we have, 

N(t) « t", P(t) » t- s 7 R 2 {t) fts t z \ 

Also, 8 — (3/vz. 
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Fig. 4 

We show the spreading of turbulent (white) sites in a laminar (black) 
background with 2 initial seeds. The horizontal axis is space and time is along 
the vertical axis. Q = 0.049, e = 0.8494. Lattice size, L = 500. Every 10 i/l 
iterate has been plotted. 

The computed values for all the exponents above defined at the 
critical values fl = 0.049, e = 0.8494 and fl = 0.073, e = 
0.4664 are listed in Table I. It is clear that the calculated ex- 
ponents agree very well with the exponents for directed per- 
colation. Thus the spatio-temporal intermittency at the points 
indicated by crosses in the phase diagram (Fig.l) belongs con- 
vincingly to the directed percolation class. 

TABLE I 

This table compares static CML exponents with the static DP exponents 



Q 


e 


z 


P/vz 


(3 


V 


rf 


0.049 


0.8494 


1.578 


0.16 


0.278 


1.1 


1.51 


0.073 


0.4664 


1.58 


0.16 


0.278 


1.1 


1.515 


DP 




1.58 


0.16 


0.28 


1.1 


1.51 



TABLE II 

This table compares dynamic CML exponents with the dynamic DP exponents 



n 


e 


V 


6 


Zs 


0.049 


0.8494 


0.313 


0.16 


1.26 


0.073 


0.4664 


0.313 


0.16 


1.26 


DP 




0.313 


0.16 


1.26 



V. Discussion 

It is interesting to note that the parameter values which cor- 
respond to spatio-temporal intermittency all lie close to the bi- 
furcation boundaries where bifurcations from synchronised so- 
lutions take place. The uppermost point where spatio-temporal 
intermittency is seen corresponds to a tangent-period doubling 
bifurcation, where two eigen-values of the stability matric cross 
+1 and —1 respectively, whereas a tangent-tangent bifurcation 
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where several of the eigen values cross the unit circle takes place 
at the other values. The synchronised solutions go to spatio- 
temporally chaotic solutions via spatio-temporal intermittency 
at these values, whereas they go to spatio-temporal chaos via 
cluster solutions in other places. Spatio-tempral intermittency 
can also be seen at other locations in the phase diagram. The 
exact nature of this behaviour is presently under investigation. 
We also hope to look at these and other aspects of dynamical 
behaviour in further detail in future work. 
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